Modeling  the  SAR  Signature  of  Nonlinear  Internal  Waves 


FINAL  REPORT 

Grant  Number:  N000 14-05- 1-0368 


Ellen  E.  Lettvin 
Applied  Physics  Laboratory 
University  of  Washington 

OBJECTIVES  AND  IMPORTANCE 

Nonlinear  Internal  Waves  are  pervasive  globally,  particularly  in  coastal  waters.  The 
currents  and  displacements  associated  with  internal  waves  influence  acoustic  propagation 
and  underwater  navigation,  as  well  as  ocean  transport  and  mixing.  Synthetic  Aperture 
Radar  (SAR)  imagery  can  reveal  the  surface  manifestations  of  internal  waves  (IWs)  in 
satellite  imagery  and  so  is  routinely  used  to  locate  and  characterize  these  features. 

Though  some  of  the  mechanisms  that  link  the  SAR  signatures,  surface  processes,  and  the 
underlying  internal  structures  have  been  understood  for  decades,  a  complete 
characterization  has  yet  to  emerge,  making  SAR  imagery  useful  only  as  a  qualitative  tool. 
The  objective  of  this  research  is  to  develop  and  validate  a  forward  model  to  predict  the 
SAR  signature  of  NLIWs  that  explicitly  includes  relevant  mechanisms  that  impact  the  sea 
surface  roughness  and  corresponding  backscattering  cross  section,  such  as  wind  speed 
and  direction,  compound  modulation  (i.e.  modulation  of  intermediate-scale  waves  by 
IWs,  which  in  turn  modulate  smaller  waves),  microscale  breaking  and  breaking  waves. 

SIGNIFICANT  RESULTS  AND  LESSONS  LEARNED 


1.  Background 


For  the  first  part  of  this  effort  we  developed  a  computationally  efficient,  web-based 
implementation  of  the  Lyzenga  and  Bennett  (1988)  (LB)  model,  which  uses  the  action 
balance  equation  to  model  the  interactions  between  surface  waves  and  currents,  such  as 
those  generated  by  the  passage  of  an  IW.  The  action  balance  equation  can  be  expressed, 
in  general  as: 
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where  N(k)  is  the  wave  action  spectrum  (see  Phillips,  1980)  at  a  given  wavenumber,  k; 
Q(k)  is  the  action  source  term  that  includes  contributions  from  wind  input,  wave 
interactions  and  dissipative  processes;  x  indexes  spatial  location  along  the  current  profile 
attributable  to  an  IW,  to  is  the  radian  frequency,  cgj  is  the  group  velocity  and  Uj  is  the 
current  due  to  the  internal  wave  in  the  direction  indexed  by  x. 
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In  LB,  the  action  balance  equation  is  governed  by  a  source  function  that  considers  energy 
input  from  the  wind  and  viscous  dissipation,  of  the  form: 


-P-(A-Ao), 
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where  /?  =  0.04|^ —  o\cos{<p  -  <pj\  is  the  wave  growth  relaxation  rate  from  Plant  (1982), 

u*  is  friction  velocity,  c  is  phase  velocity  (c=oVk),  co  is  radian  frequency  and  <j>w  is  the 
wind  direction.  The  term  Ao  is  the  equilibrium,  or  unperturbed  spectrum  for  wind¬ 
generated  surface  waves.  A  refers  to  the  perturbed  surface  waves  that  have  been 
modulated  by  an  internal  wave.  The  ratio  of  A/Ao  can  be  used  to  characterize  the 
modulation  of  surface  waves  by  currents  from  IWs.  Examples  of  A/Ao  for  waves  of 
various  scales,  ranging  from  centimeter-scale  to  tens  of  meters,  for  a  given  current  profile 
are  shown  in  Figures  1  and  2. 


For  the  LB  model,  the  equilibrium  spectrum  Ao  is  given  by: 
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A„  ( kx ,  ky)  =  N(kx,  ky)  /  D(kx,  ky) 
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where  B  is  a  constant,  ko  is  radar  wavenumber,  kx  and  ky  correspond  to  wavenumber, 
aligned  with  and  orthogonal  to  the  IW  propagation  direction. 

The  initial  implementation  of  this  model  considered  predictions  of  spectral  modulation 
for  various  choices  of  wind  speed  and  direction,  IW  propagation  speed,  direction  and 
current  profile.  These  results  were  compared  with  results  obtained  by  LB  (Lyzenga, 
personal  communication).  Scenarios  obtained  in  this  investigation  include  aligned  and 
opposing  wind  and  IW  propagation,  as  well  as  oblique  angles  of  interaction.  Our  model 
considers  both  traditional  ‘relaxation’  mechanisms  (i.e.  converging  and  diverging  surface 
currents  due  to  IWs),  as  well  as  modulation-related  mechanisms  and  associated  features 
(e.g.  breaking  and  microscale  breaking  waves). 

The  spectral  and  modulation  models  discussed  in  this  paper  are  available  on  the  www  for 
community-wide  access  (httDs://svn.anl.washington.cdu/packages/sar  iwl.  though  to  access  the 
site  requires  permission.  This  provides  a  useful  platform  for  testing,  verifying  and 
comparing  various  Equilibrium  Spectral  models,  as  well  as  various  models  that  express 
internal  wave-surface  wave  (IW-SW)  interactions.  We  have  implemented  three  spectral 
models  and  two  models  for  IW-SW  modulation;  both  are  described  in  some  detail  in  the 
literature,  so  we  provide  only  a  brief  overview  in  this  paper. 
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2.  Validation  Database 


To  assist  with  model  development  and  validation  we  assembled  a  database  that  includes 
collocated  satellite  data  (ENVISAT,  ERS,  SAR  imagery,  SPOT  imagery,  and  QuikScat 
wind  vector  maps),  ship-based  radar  measurements  (from  the  US-based  R/V  Revelle  and 
the  Taiwanese  ship  OR3),  and  moored  thermistor  chain  (S7)  measurements  acquired  in 
the  South  China  Sea  during  the  intensive  data  collection  conducted  in  the  summer  of 
2005  and  off  the  coast  of  New  Jersey  in  the  summer  of  2006.  This  database  is 
complementary  to  the  NLIWI  project  database  (nliwi.org)  and  is  one  component  of 
efforts  to  develop  a  forward  model  for  predicting  the  sea  surface  patterns  attributable  to 
IWs  and  their  corresponding  SAR  signatures. 

IW  events  were  found  in  all  observational  data  including  those  acquired  by  ships,  moored 
platforms  and  satellites  including  ENVISAT,  ERS-2,  and  SPOT.  Though  IW  events 
could  be  readily  identified  in  individual  data  sets,  there  was  a  surprisingly  limited  number 
of  cases  where  there  were  simultaneous  ship-based  and  satellite-based  or  moored 
observations  of  IWs.  A  likely  reason  for  the  lack  of  collocations  is  that  ship-based  radars 
have  a  limited  range  (0(10  km))  relative  to  the  scale  of  satellite  imagery,  so  though 
considerable  SAR  imagery  containing  NLIWs  were  acquired  during  both  experiments, 
the  number  of  collocations  is  extremely  limited  due  to  constraints  including  time  and 
location  of  satellite  overpasses,  direction  of  ship  heading  and  radar  look  direction.  Figure 
3  provides  an  example  of  a  collocated  satellite-  and  ship-based  observation  of  IWs;  the 
ship-based  radar  data  were  acquired  from  the  Taiwanese  Research  Vessel,  the  OR3. 


IW  Current  ProlUe,  Constant  In  y  Direction;  Wind  Speed  =4  m/s.  Wind  Dir  =  180  deg  CW  from  x-Axis 


IW  Modulation  of  Surface  Waves  Propagating  180  Deg  from  x-Axis;  IW  phase  spd  =  2  m/s  in  +x  Direction 


Figure  1 .  Modulation  of  equilibrium  spectral  density  by  currents  due  to  NLIWs.  Upper  panel  shows  the 
surface  currents  in  the  direction  of  propagation  of  the  NLIW;  wind  speed  is  4ms'1,  and  the  NL1W 
propagation  is  directly  opposing  the  wind.  Lower  panel  shows  the  modulation  (A/Ao)  experienced  by 
surface  waves  of  1  cm  (*)  and  10  m  (o)  due  to  interaction  with  the  NLIW;  1cm  (Bragg-scale)  waves 
experience  negligible  modulation,  whereas  10m  waves  experience  significant  modulation. 
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IW  Current  Profile.  Constant  m  y  Direction;  Wind  Speed  =4  m/s.  Wind  Dir  =  180  deg  CW  from  x-Axis 


IW  Modulation  of  Surface  Waves  Propagating  180  Deg  from  x-Axis;  IW  phase  spd  *  2  m/s  in  ♦x  Direction 


Figure  2.  Modulation  experienced  by  waves  with  wavelengths  between  the  two  extremes  presented  in 
Figure  1 .  Upper  panel  shows  the  surface  current  distribution  for  same  conditions  as  in  Figure  1 .  Lower 
panel  shows  the  modulation  (A/Ao)  experienced  by  surface  waves  of  3  cm  (*),  10  cm  (o),  30  cm  (square), 
lm  (triangle)  and  4  m  (diamond)due  to  interaction  with  the  NLIW.  Again  the  Bragg-scale  waves  (<  10  cm) 
experience  negligible  modulation,  whereas  the  waves  of  0.3,  1  and  4m  experience  significant  and 
progressively  increasing  modulation. 

3.  Implementation 

Our  first  step  was  to  implement  and  then  validate  the  LB  model  via  qualitative  and 
quantitative  comparisons,  which  we  accomplished  successfully.  Following  preliminary 
implementation  and  validation,  undesirable  features  were  noted  in  the  LB  spectral  model, 
particularly  a  singularity  found  at  +/-  n!2  relative  to  the  wind  direction.  To  avoid  this 
feature  we  implemented  an  alternative  version  of  this  scheme  that  used  the  spectral  model 
of  Donelan  and  Pierson  (1987)  (hereafter  referred  to  as  DP)  in  the  action  modulation 
framework  developed  by  LB.  This  model  was  selected  for  the  following  reasons:  it  was 
derived  based  on  physical  principles  and  measurements  and  has  been  used  in  many 
studies  in  radar  remote  sensing  research.  Also,  at  low  wavenumbers,  the  DP  model  is 
identical  to  spectral  model  of  Kudryavtsev  et  al.(  1 999),  which  is  used  in  a  subsequent 
part  of  this  effort.  The  DP  spectrum  employs  the  same  source  terms  and  action 
modulation  framework  for  computing  SW-IW  interactions  as  the  LB  model,  but  is  more 
physically  realistic. 

The  inputs  to  the  LB  model,  for  implementation  of  either  the  LB  or  DP  spectral  model 
includes:  wind  speed  and  direction,  IW  speed  and  direction,  maximum  surface  current, 
current  distribution,  and  current  extent.  Given  these  inputs,  the  LB  model  obtains  the 
two-dimensional  wavenumber  spectra  (in  terms  of  energy,  action  or  saturation,  where 
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Figure  3.  This  figure  shows  collocated,  satellite-based  SAR  and  ship-based  radar  measurements  of 

an  internal  wave  in  the  South  China  Sea.  This  C-band,  vertically  polarized  ENVISAT  image  was  acquired 
April  28,  2005;  wind  speed  was  0-2  ms'1  and  from  the  SE,  as  estimated  from  a  contemporaneous 
QuikScat  overpass.  The  location  of  the  OR3  is  indicated  in  red  in  the  upper  left  comer  of  the  image.  The 
three  panels  on  the  right  of  the  figure  show  (top)  an  X-band  ship-based  radar  measurement  obtained  from 
the  OR3;  (middle)  an  intensity  plot  taken  in  the  direction  of  the  IW  (indicated  by  a  red  line  in  the  top  panel) 
that  clearly  shows  the  presence  of  an  IW;  (bottom)  an  intensity  plot  taken  from  the  ENVISAT  image  along 
the  direction  indicated  by  the  red  line  originating  at  the  location  of  the  OR3.  The  red  shaded  area  shows  the 
portion  of  the  satellite-based  intensity  plot  overlapping  with  the  intensity  plot  for  the  ship-based  radar. 
Though  the  resolution  of  the  ENVISAT  image  is,  as  expected,  much  coarser  than  the  ship-based  radar,  the 
same  IW  can  still  be  clearly  observed  in  both  plots. 


saturation,  B,  is  given  by  B(kx,ky)  =  k4S(kx,ky);  k  is  wavenumber  and  S(kx,ky)  is  the 
wavenumber  spectrum,  and  is  a  function  of  x-  and  y-wavenumber,  as  well  as  position 
along  the  internal  wave,  x,  reflecting  the  modulation  of  the  surface  wave  action  with  the 
passage  of  the  IW. 
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Don  el  an  Equilibrium  Saturation  Spectrum 
Wind  Speed  =  7  m/s 


Figure  4  This  shows  the  equilibrium  spectral  model  of  Donelan  and  Pierson  (1987)  for  a  wind  of 

7ms'1  blowing  from  right  to  left;  the  color  bar  shows  spectral  levels  in  dB. 

An  example  of  the  Donelan  and  Pierson  Equilibrium  spectrum  model  is  shown  in  Figure 
4.  Here,  the  wind  speed  is  7ms'1,  wind  direction  is  from  the  left  to  the  right  of  the  figure 
(from  7t),  and  the  equilibrium  spectrum  is  depicted  in  terms  of  saturation. 

In  addition  to  the  LB  scheme,  we  implemented  an  approach  proposed  by  Kudryavtsev  et 
al.  (2003,  2005)  for  predicting  IW-SW  modulation.  This  approach  uses  the  innovative 
spectral  model  of  Kudryavtsev  et  al.  (1999)  and  includes  an  action  source  term  associated 
with  the  generation  of  parasitic  capillary  waves  produced  from  the  IW-SW  interaction. 

An  illustration  of  this  spectral  model  is  provided  in  Figure  5,  for  the  same  conditions 
(wind  speed  and  direction)  used  to  illustrate  the  DP  model  (Figure  4). 

The  action  source  terms  used  by  Kudryavtsev  (2003)  are  given  in  Equation  5. 


Q(k)  =  0)ik~5 


pv{k)B{k)-B{k) 


Bjk) 

a  ) 


+  V(*) 


(5) 


The  first  term,  pv(k)  =  P(k)  -  4vk2/co  is  ‘effective  growth  rate’,  defined  as  wind  growth 
rate  less  viscous  dissipation,  where  P  is  given  by 

P(k)  =  Cp(u*/c)2  cos2  <p  (6) 
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(p  is  the  angle  between  wind  and  wave  direction,  and  Cp  is  given  by  Stewart  ( 1 974). 

The  second  term  of  (5)  expresses  nonlinear  wave-wave  interaction;  a  and  n  are  functions 
of  k/kg,  where  ky  =(g/y)u 2  =  k  of  minimum  phase  velocity,  and  kg=  ky2/k; 

The  third  term  in  expression  (5)  Ipc  expresses  the  contribution  of  parasitic  capillary  waves 
to  action  modulation  and  is  given  by 

Ipc(k)  =  Pv(kg)B(kg)<|>(k/kg),  (7) 

where  <J>(k/kg)  is  a  filter  function  which  restricts  the  action  in  k-space  of  source  term  Ipc(k) 
to  those  gravity  waves  that  generate  parasitic  capillaries.  Figure  6  illustrates  how  1^ 
works.  The  upper  two  panels  show  the  wavenumber  dependence  of  the  terms  a  and  n  in 
the  second  action  source  term  (in  Equation  6).  The  bottom  panel  illustrates  the 
wavenumber  dependence  of  Ipc,  which  is  essentially  a  filter  function.  Ipc  will  be  nonzero 
when  conditions  are  favorable  for  the  generation  of  parasitic  capillary  waves  by  gravity 
waves. 


Kudryavtsev  Equilibrium  Saturation  Spectrum 
Wind  Speed  =  7  m/s 


Figure  5  This  shows  the  equilibrium  spectral  model  of  Kudryavtsev  (1999)  for  a  wind  of  7ms'1 

blowing  from  right  to  left;  the  color  bar  shows  spectral  levels  in  dB. 


Variation  of  the  Ipc  source  term  with  wavenumber  is  shown  in  Figure  7,  which  also 
illustrates  the  relative  contribution  of  Ipc  to  the  Kudryavtsev  spectrum.  The  panel  on  the 
left  shows  the  two  dimensional  wavenumber  spectral  model  of  Kudryavtsev  et  al.  (1999) 
for  a  wind  of  7ms- 1 ,  in  the  n  direction  (from  right  to  left);  the  panel  on  the  right 
illustrates  the  magnitude  of  the  source  term  for  parasitic  capillary  waves,  Ipc,  as  a  function 
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Figure  6  .  Illustration  of  the  wavenumber  dependence  of  the  exponent  a  (upper  panel),  and  denominator  n 
(middle  panel),  in  the  Kudryavtsev  (2003)  model’s  action  source  term  for  wave-wave  interactions.  The 
bottom  panel  shows  the  filter  corresponding  to  the  presence  of  parasitic  capillary  waves.  The  red  and  blue 
lines  are  included  to  show  the  relationship  between  gravity  waves  and  corresponding  parasitic  capillary 
waves. 


Kudryavtsev  Equilibrium  Saturation  Spectrum 
Wind  Speed  =  7  m/s 


Figure  7  .  Kudryavtsev  (2003)  equilibrium  spectral  model  (saturation  form)  and  corresponding 

action  term  associated  with  parasitic  capillary  waves,  1^.  Color  bar  shows  relative  power  levels  in  dB. 
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of  wavenumber  and  direction.  Note  that  the  range  of  wavenumbers  illustrated  in  Figure  7 
is  expanded  relative  to  what  is  depicted  in  Figure  5;  the  color  bar  (in  dB)  is  also  scaled  to 
match  the  scaling  of  the  source  term  associated  with  parasitic  capillary  waves,  I,*. 

We  compare  the  DP  and  Kudryavtsev  (1999)  equilibrium  spectral  models  (i.e.  not 
including  SW-IW  modulation  effects)  in  Figure  8.  Though  not  included  in  the  analyses 
described  here,  it  is  worth  noting  that  the  2-dimensional  wavenumber  spectrum  serves  as 
the  inputs  to  electromagnetic  scattering  models  used  to  predict  the  SAR  signature  of  the 
sea  surface.  We  intend  to  implement  the  stochastic,  multiscale  model  of  Plant  (2004), 
and  will  compare  model-based  predictions  with  satellite-  and  ship-based  radar 
measurements  as  part  of  model  validation.  Recall  that  the  two  approaches  (LB  action 
computation  with  DP  spectral  model,  and  Kudryavtsev  modulation  and  spectral  model) 
use  different  action  source  terms  and  different  computational  schemes. 


Donelan  Equilibrium  Saturation  Spectrum 
Wind  Speed  =  7  m/s 


kx 


Kudryavtsev  Equilibrium  Saturation  Spectrum 
Wind  Speed  =  7  m/s 


-2000  0  2000 
kx 


Figure  8  2-Dimensional  wavenumber  spectral  model  of  Donelan  et  al.  (1985),  left  panel,  and 

Kudryavtsev  et  al.  ( 1 999),  right  panel,  for  wind  speed  7ms'1;  wind  direction  =  7t.  Color  bar  shows  relative 
power  in  dB.  Note  color  bar  is  different  in  Figures  7  and  8. 

Figures  9  and  10  enable  us  to  compare  the  modulated  spectrum  at  various  locations  along 
a  current  profile  produced  by  a  hypothetical  IW  for  both  the  DP/LB  model  and  the 
Kudryavtsev  model.  In  these  examples,  we  assume  the  current  profile  has  a  sech2  shape, 
extends  over  3  km  (from  -1 500  m  to  +  1 500  m),  and  attains  a  maximum  of  0.7  ms'1,  as 
illustrated  in  the  bottom-most  panel  of  each  figure.  The  indices  A  -  D  assigned  to  the 
upper  panels  of  the  figure  correspond  to  location  along  the  current  profile  shown  in  the 
lower  panel.  The  spectrum  at  location  D  will  be  similar,  but  not  identical  to  the 
equilibrium  (unmodulated)  spectrum,  Ao. 
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Figure  9.  Upper  four  panels  illustrate  the  impact  of  IW-SW  modulation  on  surface  waves 
distributed  according  to  the  Donelan  and  Pierson  (1987)  spectral  model.  Model  results 
are  provided  corresponding  to  surface  waves  between  62  m  2.5  m  wavelength.  We  focus 
here  on  intermediate-scale  waves,  as  the  LB  model  assumes  (as  shown  in  Figures  1  and 
2)  that  IWs  directly  modulate  centimeter-scale  (Bragg)  waves.  The  color  bar  shows  the 
saturation  spectrum  (in  dB)  corresponding  to  the  four  locations  indexed  by  letters  A-D. 
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-2000  -1500  -1000  -500  0  500  1000  1500  2000 


Figure  10.  As  in  Figure  9,  the  upper  four  panels  illustrate  the  impact  of  IW-SW 
modulation  on  surface  waves  distributed  according  to  the  Kudryavtsev  ( 1 999)  spectral 
model.  Model  results  are  provided  corresponding  to  surface  waves  between  62  m  2.5  m 
wavelength.  We  focus  here  on  intermediate-scale  waves,  as  the  Kudryavtsev  model  also 
assumes  that  IWs  directly  modulate  centimeter-scale  (Bragg)  waves.  The  color  bar 
shows  the  saturation  spectrum  (in  dB)  corresponding  to  the  four  locations  indexed  by 
letters  A-D. 
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During  the  field  experiment  held  off  the  coast  of  New  Jersey  in  the  summer  of  2006, 
surface  current  velocity  estimates  were  obtained  from  the  R/V  Oceanus  (J.  Mourn, 
personal  communication)  that  revealed  the  current  structure  shown  in  Figure  11.  An 
analytical  fit  (D.  Lyzenga,  personal  communication)  was  obtained  for  these  data,  which 
corresponds  to  the  red  line  superimposed  on  the  blue  data  plot.  It  is  believed  that  the 
sech2-type  current  profile  shown  in  Figures  9  and  1 0  corresponds  to  I W  characteristics 
typically  observed  in  the  South  China  Sea;  current  profiles  such  as  shown  in  Figure  1 1 
correspond  to  what  is  frequently  observed  off  the  coast  of  New  Jersey.  In  Figures  12  and 
1 3  we  look  at  the  I W-S W  modulation  in  different  wavenumber  regimes  produced  by  a 
current  profile  similar  to  that  shown  in  Figure  1 1 . 
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Figure  1 1  Surface  current  estimates  for  a  train  of  internal  waves  observed  off  the  coast  of 
New  Jersey  in  summer  2006,  extrapolated  to  the  surface  from  ADCP  measurements  on 
the  R/V  Oceanus  (courtesy  of  J.  Mourn  for  observational  data,  and  D.  Lyzenga  for  the 
extrapolation).  Blue  lines  indicate  observed  data;  heavy  red  lines  correspond  to  a 
polynomial  fit  to  the  data;  coefficients  noted  on  the  lower  edge  of  the  plot.  These  current 
estimates  are  used  in  Figures  12  and  13  to  modulate  the  surface  wave  spectra. 
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Figure  12  Here  we  can  examine  IW-SW  modulation  of  the  Kudryavtsev  spectrum  for 
wavenumbers  between  0.1  and  2.5  (wavelengths  between  62.5  m  and  2.5  m).  Panels  A  - 
D  show  modulated  SW  spectrum  at  various  locations  along  an  IW  current  profile  (shown 
in  bottom  panel);  letters  designate  locations  where  spectra  are  obtained.  Color  bar  shows 
spectral  saturation  levels  in  dB. 


13 


_]L  I  I  I  I  I  I  I  I 

-1200  -1000  -800  -600  -400  -200  0  200  400 


Figure  13  Here  we  consider  IW-SW  modulation  of  the  Kudryavtsev  spectrum  for  a 
much  broader  range  of  wavenumbers  between  0.1  and  500  (wavelengths  between  62.5  m 
and  0.01  m).  Panels  A  -  D  show  modulated  SW  spectrum  at  the  same  locations  along  an 
1W  current  profile  as  were  considered  for  Figure  12;  letters  designate  locations  where 
spectra  were  obtained;  again  locations  are  indicated  in  bottom  panel.  Color  bar  shows 
spectral  saturation  levels  in  dB. 


14 


SUMMARY 


We  have  successfully  applied  the  spectral  models  of  Lyzenga  and  Bennett  (1988), 
Donelan  and  Pierson  (1987)  and  Kudryavtsev  et  al.  (1999, 2003),  as  well  as  the  action 
modulation  frameworks  of  Lyzenga  and  Bennett  and  Kudryavtsev  (2003).  We  used  the 
various  techniques  to  study  the  modulation  of  surface  waves  by  internal  waves.  As 
Figures  9-13  demonstrate,  we  were  able  to  show  that  there  is  considerable  IW-SW 
modulation,  but  none  of  the  approaches  implemented  thus  far  are  completely  satisfactory. 
For  example,  the  model  of  Kudryavtsev  et  al.  assumes  that  the  longer  waves  produce 
parasitic  capillary  waves,  but  small-scale  breakers  are  considered  to  be  sources  of  short 
waves  in  the  action  conservation  equations.  But  these  equations  are  usually  solved 
assuming  that  the  short  waves  move  at  their  intrinsic  phase  speeds,  an  assumption  that  is 
not  correct  for  short  waves  that  are  being  advected  by  longer  ones,  let  alone  for  those  that 
are  forced  by  the  longer  waves.  Furthermore,  models  for  predicting  backscatter  require 
us  to  know  where  on  longer  waves  the  breakers  and  parasitic  capillaries  are  located,  and 
this  cannot  be  determined  with  the  present  approach. 


NEXT  STEPS 

While  the  work  described  in  this  paper  has  made  considerable  progress  towards 
understanding  IW-SW  modulation,  and  towards  predicting  the  surface  manifestations  of 
this  phenomenon,  there  is  clearly  considerable  work  still  to  be  done.  Additional  physical 
effects  need  to  be  incorporated  into  the  spectral  models  and  action  modulation  framework 
models.  Plant  has  proposed  to  follow  on  the  research  described  in  this  report  with  a  new 
approach  where  parasitic  capillaries  and  breaking  wave  patches  are  considered  distortions 
of  longer  waves  in  the  scattering  model,  rather  than  being  included  in  the  action 
conservation  equations.  Also,  a  more  recent  paper  by  Kudryavtsev  et  al.  (2005)  considers 
action  modulation  to  arise  from  wind  forcing,  viscous  effects,  dissipation  via  wave 
breaking,  wave-wave  interactions,  generation  of  free  waves  by  breaking  waves,  and 
generation  of  bound  (parasitic)  capillary  waves  by  breaking  waves.  While  this  approach 
looks  promising,  it  does  not  address  the  concerns  raised  above  regarding  short  wave 
advection  by  longer  waves. 

Additional  next  steps  we  hope  to  validate  the  models  described  in  this  report  are  to 
compare  predictions  of  IW-SW  modulation  with  measurements  from  ASIS  buoys  (Hans 
Graber,  personal  communication)  deployed  during  field  experiments  in  regions  where 
significant  IW  activity  is  anticipated,  e.g.  South  China  Sea  or  off  the  coast  of  New  Jersey. 
In  addition  to  comparisons  of  direct  measurements,  we  would  like  to  compare  predictions 
of  a0  (radar  backscatter)  based  on  with  ship-  tower-  or  satellite-based  radar  measurements 
for  a  variety  of  current  configurations.  Final  improvements  we  hope  to  implement 
include  prediction  of  modulation  from  currents  that  vary  both  in  x  and  y;  interactions 
involving  SW,  IW  and  swell,  and  implementation  of  a  faster,  more  efficient  version  of 
the  code. 
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